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ABSTRACT 



Context. The numerous results obtained with asteroseismology thanks to space missions such as CoRoT and Kepler are providing a 
new insight on stellar evolution. After five years of observations, CoRoT is going on providing high-quality data. We present here the 
analysis of the double star HD 169392 complemented by ground-based spectroscopic observations. 

Aims. This work aims at characterizing the fundamental parameters of the two stars, their chemical composition, the acoustic-mode 
global parameters including their individual frequencies, and their dynamics. 

Methods. We have analysed HARPS observations of the two stars to retrieve their chemical compositions. Several methods have been 
used and compared to measure the global properties of acoustic modes and their individual frequencies from the photometric data of 
CoRoT. 

Results. The new spectroscopic observations and archival astrometric values suggest that HD 169392 is a wide binary system weakly 
bounded. We have obtained the spectroscopic parameters for both components, suggesting the origin from the same cloud. However, 
only the mode signature of HD 169392 A has been measured within the CoRoT data. The signal-to-noise ratio of the modes in 
HD 169392B is too low to allow any confident detection. We were able to extract mode parameters of modes for {=0, 1, 2, and 3. 
The study of the splittings and inclination angle gives two possible solutions with splittings and inclination angles of 0.4-1.0yuHz and 
20-40° for one case and 0.2-0.5 yuHz and 55-86° for the other case. The modeling of this star with the Asteroseismic Modeling Portal 
led to a mass of 1.15±0.01M o , a radius of 1. 88 + 0.02 Rg, and an age of 4.33 ±0.12Gyr, where the uncertainties are the internal 
ones. 

Key words. Asteroseismology - Methods: data analysis - Stars: oscillations - Stars: individual: HD 169392 



1. Introduction 

The convective motions in the external layers of solar-like oscil- 
lating stars excite acoustic sound waves, which become trapped 
in the stellar interiors (e.g., Goldreich & Keeley]|1977 Samadi 



2011 for a detailed review). The precise frequencies of these 
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* The CoRoT space mission, launched on December 27th 2006, 
has been developed and is operated by CNES, with the contribution 
of Austria, Belgium, Brazil, ESA (RSSD and Science Programme), 
Germany and Spain. 

** This work is based on ground-based observations made with the 
ESO 3.6m-telescope at La Silla Observatory under the ESO Large 
Programme LP185-D.0056. 



pressure (p) driven waves depend on the properties of the 
medium in which they propagate. Thus, stellar seismology al- 
lows us to infer the properties of the Sun and stellar interiors by 
studying and characterizing these p modes (e.g., |Gough et al 



1996; Christensen-Dalsgaard 2002). However, our knowledge 
of the properties of the stellar interiors depends on our ability 
to correctly measure the oscillations modes (e.g. Appourchaux 
|et al |1998[|Appourchaux|2TJTT] l. 

In recent years the French-led satellite, Convection, Rotation 
and planetary T ransits (CoRoT,|Baglin et al.|2006|), an d NASA's 
Kepler mission ( Koch et al. 2010||Borucki et aL 2010 1 have been 
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providing high-quality, long-term seismic data of solar-like stars 
CoRoT has studied some interesting F- and G-type stars (e.g. 
Appourchaux et al.|2008| |Barban et al.|2009l |Garcfa et al.|2 009 



Deheuvels et al.|2010|l in cluding some hosting planets (Gaulme 



et al.||2010t |Ballot et al-l|201 lb[ |Howell et al.|[20T2] > and oth- 
ers in binary systems ( Mathur et al.|2010a| l. Kepler observations 
have enabled ensemble asteroseismology of hundreds of solar- 
like stars (e.g. Chaplin et al.|201 fb| ) as well as precise studies of 
long time series with more than 8 months of nearly continuous 



data (ICampante et aL]|201 1 Mathur et al.|2011a Appourchaux 



et al.|2012b| l, and pulsating stars in binary systems (e.g. Hekker 
et al.[2010bf White et al. in prep.) and in clusters (e.g. Basu et al. 
2011||Stello et al.||2011b| l. The fundamental properties of stars 
can be inferred by means of asteroseismic quantities either by 
using scaling relations (e.g. Kjeldsen & Bedding||l995| 2011 



Samadi et al.|2007b|l, or by using stellar models (e.g.|Piau et aL 
2009}|Metcalfe et al.|2010||20"lfj|Creevey et al.|2012[|Escobar 
et al.|2012||Mathur et al.|2012||Mazumdar et al.|2012) . This al- 
lows us to have some information on the structure of the stars 
but also to test the models and start improving the physics in- 
cluded in stellar evolution codes (e.g. Christensen-Dalsgaard & 
Thompson 201 1 and references there in). 

In the present paper we report the study of the double 
star HD 169392, using high-resolution spectrophotometry by 
HARPS combined with seismic analysis from CoRoT photom- 
etry. However in this latter case, only HD 169392 A was ana- 
lyzed because the HD 169392B was too faint. The two com- 
ponents of HD 169392= WDS 18247-0636 have a separation 
of 5.881+0.006", a position angle of 195.08±0.09 °, and a 
magnitude difference A V- 1.507 ±0.002 at the epoch 2005.48 



(Sinachopoulos et al. 2007). The two components are reported in 
SIMBAD^jas having similar proper motions. HD 169392A=HIP 
90239 is the main component of the pair, with spectral type 
GO IV and a parallax n = 13.82 ± 1.20 mas (van Leeuwen 
2007). From the Geneva-Copenhagen survey (Nordstrom et al. 
2"004"j |Holmberg eTaLl[2007] l this primary star has T eS = 5942 
K, [Fe/H] = -0.03, m v = 7.50 mag, and vsin; = 3 kms" 1 . The 
second component of the system, HD 169392B ( Duncan|1984| l 
is a GO V - G2 IV star with m v = 8.98 mag. 

New spectroscopic observations performed by HARPS 
( |Mayor et al.| |2003 ) and detailed in Section [2] have allowed a 
precise characterization of the stellar properties of the two stars 
in the system. In Section 3, we describe the methodology used 
to prepare the original time series from CoRoT, to fit the mode 
oscillation power spectrum, and to derive a single frequency set 
from results of different fitters. We derive the global properties 
of the modes of both stars (Section 4) while in Section 5 we anal- 
yse the surface rotation and the background of HD 169392A. In 
Section 6, we present the peak-bagging methods used to obtain 
the frequencies of the individual p modes of HD 169392 A and 
the results are discussed in Section 7. With the spectroscopic 
constraints derived and the frequencies of modes, we modeled 
the star and retrieve the fundamental parameters of HD 169392 A 
in Section 8. Finally, the main results are summarized in Section 
9. 



2. Fundamental stellar parameters 

High-quality spectra of HD 169392 were obtained in the frame- 
work of the CoRoT Ground-Based Observations Working Group 
(e.g. |Catala et aL|2006t |Uytterhoeven et al.||2008| |Poretti et aL 



Table 1. Table with spectroscopic parameters for both A and B 
components of the system derived from our HARPS observa- 
tions and the seismic logg. 





A 


B 


Units 




5985 ± 60 


5885 ± 70 


K 




3.96 ± 0.07 


3.75 ±0.07 


cgs 


[Fe/H] 


-0.04 ±0.10 


-0.10 ± 0.10 


dex 


^micro 


1.41 ±0.07 


1.04 ±0.07 


km s -1 


v sin i 


1.0 ±0.6 


3.0 ± 0.7 


km s -1 


"macro 


2.8 ± 0.4 


2.5 ±0.4 


km s -1 (fixed) 



2012). We obtained high signal-to-noise ratio (S/N) spectra of 
each component using HARPS in the high-resolution mode 
(R=l 14,000). The spectrum of HD 169392A taken on the 25- 
26 June 201 1 night shows a radial velocity (RV) of -68.1606 ± 
0.0004 km s 1 , that of HD169392B taken on the 24-25 June 201 1 
shows RV=-70.1388 + 0.0004 kms -1 . To have common proper 
motions and compatible radial velocities is a sufficient condition 
to consider a double star as a physical pair (e.g. |Struve et al 



[T955l >. In the case of HD 169392, the separation of 0.0588 "at a 
distance of 723 pc corresponds to 4250 AU and then the system 
should be weakly gravitationally bound. 

The HARPS spectra have both an excellent S/N, i.e., 
S/N=267 for the A component, and S/N=234 for the B on e. They 
were reduced using a MIDAS semi-automatic pipeline |Rainer| 
(2003). Their high quality allowed a careful and detailed analy- 
sis by means of the semi-automatic pipeline VWA ( Bruntt et al.| 
2010). More than 400 spectral lines have been analyzed in each 



star. The match of the synthetic line profile to the observations 
was inspected by eye and several lines were rejected either due to 
problems with the placement of the continuum or due to strong 
blends in the line wings. Atmospheric models were determined 



by interpolating on a grid of MARCS models dGustafsson et al. 



2008) and line data were extracted from VALD (IKupka et al. 



1 This research has made use of the SIMBAD database, operated at 
CDS, Strasbourg, France. 



1999). The oscillator strengths (log gf) were adjusted relative to 
a solar spectrum as described by Bruntt et al. (2012). 

The classical analysis method consists of treating the atmo- 
spheric parameters of the effective temperature, T e g, the surface 
gravity, logg, and microturbulence, v macro as free parameters. 
The optimum solution is then selected to minimize the corre- 
lation between Fe i abundance, equivalent width, and excitation 
potential. Furthermore, we required that the mean abundance de- 
termined from neutral and ionized Fe lines agree. The final pa- 
rameters are listed in Table[T] For these parameters we calculated 
the abundances of 13 different species. We have done the anal- 
ysis with and without the seismic log g, which is available only 
for the A component (see Section 4.1). When the seismic value 
is used, it affects the log g value by +0.2 dex and T e g by +100 K. 
We also note that [Fe/H] is increased by about 0.05 dex. 

We listed the abundances relative to the Sun (AA) along 
with the number of spectral lines used in Table [2] and plotted 
them in Fig. [T] The lithium abundance has been determined 
for both components of the system: [L1/H]a = +1.42 + 0.10 
and [L i/H] B = +0.23 + 0.12 (see details in |Bruntt et IT] 
(2012)). According to the temperatures of both stars, the mea- 
sured lithium abundances are consistent with the sample ana- 
lyzed for example by Israelia net al.| ( [2004] l. 

We have adjusted the projected splittings, v sin i value by fit- 
ting the observations to synthetic line profiles for several isolated 
lines. During these fits the macroturbulence was fixed using the 
calibration for solar-type stars from Bruntt et al. (2010 1. We note 
that the metallicities of the two components are quite similar. 
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Table 2. Abundances relative to the Sun (AA) and the number of 
lines (AO used in the spectral analysis of the two components (A 
and B) of the binary system HD 169392. The abundances have 
an uncertainty of 0.08 dex for all the elements excepting for the 
lithium, which are 0.10 and 0.12 for A and B components re- 
spectively. 
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Fig. I. Abundances for 14 elements relative to the Sun for the A 
(top) and B (components) of the system. Circles and box sym- 
bols are used for the mean abundance from neutral and singly 
ionized lines, respectively. The dark grey horizontal bar marks 
the mean metallicity within 1-cr uncertainty range. The lithium 
is not plotted because it has a much higher value than the other 
elements in the A component. 



2.1. Deriving v max from spectroscopic parameters 

The proportionality between the frequency at maximum power, 
Vma x, and the acoustic cu toff frequency w as originally suggested 
by Brown et al.| ( |1991| l, developed by Kjeldsen & Bedding 
19951, and justified theoretically by Belkacem et al. (201 1 



According to theory, the large frequency spacing, Av, is propor- 
tional to the square root of the mean density of the star. The scal- 
ing relations - scaled to solar values- have also been extensively 
tested with observational data (e.g. |Huber et al.||201l| l. The re- 
lations allow us to obtain approximate values for the so-called 
mean "large spacing", (Ay), and the frequency of maximum am- 
plitude v max . 



Combining both equations, we are able to deduce v max as 
a function of the effective temperature and the log g derived in 
Table [T] for both components of the binary system. Taking into 
account the uncertainties in the observations we have in the spec- 
troscopic parameters, we obtain a range of v max for HD 169392A 
of [524, 730] yuHz, while for HD 169392B the range is [2202, 
3077] juHz. 



3. CoRoT photometric light curve 

CoRoT has observed HD 169392 (CoRoT id = 9161) continu- 
ously during 91.2 days starting on 2009 April 1 till 2009 July 2. 
These observations correspond to the third long run in the center 
direction of the galaxy (LRc03). In the present study, we have 
used the so called helreg level 2 (N2) datasets (Samadi et al. 
2007a), i.e., time series prepared by the CoRoT Data Center 



(CDC) that are regularly cadenced at 32s in the heliocentric 
frame. 

During the crossing of the South Atlantic Anomaly (SAA), 
the CoRoT measurements are perturbed (Auver gne et al.|2 009) 
and the power spectrum shows a sequence of spikes at fre- 
quencies of n x 161.7 /iHz where n is an integer. Moreover, 
other spurious peaks appear at multiples of daily harmonics at 
(n x 1 6 1 .7) ± (m x 1 1 .57) //Hz, where m is an integer. Therefore, 
some teams in the collaboration used different interpolation tech- 
niques to fill the perturbed data. One method used an "inpaint- 
ing" algorithm -a Multi-Scale Discrete Cosine Transform (Sato 
et al.|2 010)- because we obtained good results for other CoRoT 



targets (see for example HD 170987, Mat hur et aLpO lOa). The 
overall duty cycle before the interpolation was 83.4%. Another 
method was to interpolate the gaps produced by SAA with 
parabola fitting the points around the gap, in a way similar to 
the usual interpolation performed in N2 data for missing points 
(see |Ballot et al.|201 lb] for details). 

The measured light curve has been converted into a relative 
flux (ppm) by correcting first for any discontinuity in the flux 
and then removing a sixth order polynomial fit to take into ac- 
count the aging of the instrument (Auverg ne et al.|[20 09). The 
resultant flux is plotted in Fig. [2 Two jumps -of possible in- 
strumental origin- are visible at the dates 43.3 and 80.05 days. 
Although they have no significant influence in the power spectral 
density (PSD), they have an important impact on the analysis of 
the surface rotation (see Section [5T2] >. To compute the PSD, we 
used a standard fast Fourier transform algorithm and normalized 
it as the so-called one-sided power spectral density (Press et al. 
1992). Finally, we calibrated the spectrum to comply with the 
Parseval's theorem. 



4. Global seismic parameters of the binary system 

4.1. HD 169392 A 

Three different teams, A2Z, COR and OCT, analyzed the global 
seismic properties of HD 169392 A using different method- 
ologies. Explanation of each method can be fou nd in |Mathur 
et al.|(|2010b[ ), |Mosser & Appourchaux| ( |2009] l and |Hekker etal. 
(2010ab, respectively. Each team estimated the mean large fre- 
quency spacing, (Av), the frequency of maximum power in the p- 
mode bump, v max , and the maximum amplitude per radial mode, 
converted to the bolometric one following the method described 
in|Michel et al.|(|2009|l, Aboi,i=o- The results from A2Z are sum- 



marized in Table [3] Briefly, the A2Z and OCT methods analyse 
the power spectrum of the power spectrum to measure the mean 
large separation, while the COR method computes the envelope 
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a. 
a, 

x 

3 




40 60 
Time (days) 

Fig. 2. N2-helreg relative flux (in grey) after inpainting the bad 
points and those taken during the SAA crossing. The aging has 
been removed by a sixth order polynomial fit. The black curve 
corresponds to a 1 hour rebin. 

Table 3. Global seismic parameters computed using the A2Z 
pipeline. 



Quantity 


Value 


Units 


<Av> 


56.32 ± 1.17 


//Hz 


^max 


1030 ± 54 


//Hz 


■Abol,l=0 


3.61+0.35 


ppm 



autocorrelation function. All teams show good agreement on the 
values of the seismic parameters when we take into account the 
different frequency ranges used for the calculations (a full com- 
parison of them can be found in Hekker et al. ( 201 l| l and |Verner 
|etaL| ( j20TT) ). 

The values of v max and Ab i,i=o were obtained by smoothing 
the background-subtracted PSD with a sliding window of width 
2(Av) as described by Kjeldsen et al. ( 2008a| ). Then, to derive 
■Aboi,i=0 we followed the procedure to correct for the CoRoT in- 
strumental response as described in Michel et al. ( 2008 1, while 



Vmax corresponds to the frequency of the maximum power of a 
fitted Gaussian. 



Using scaling relations from solar values (e.g. Kjeldsen 
& Bedding 1995] > and the effective temperature derived by 
our HARPS observations (see Table [TJ, we can infer some 
seismic stellar parameters. Therefore we obtain a mass of 
1.34 + 0.26 M Q , a radius of 1.97 ±0. 19 R© and a logg of 3.96 
+ 0.013 dex. 



nent by 80 %, and by assuming a radius of ~ 1.03 R e and a mass 
of ~ 1 M from the scaling relations, we obtain a probability of 
2 % to detect the modes. This agrees with the non-detection with 
the CoRoT data. 

We focus the remainder of the paper on the A component of 
the system. 



5. Background and surface Rotation of HD 169392A 

5.1. Stellar background 

From the measured flux we compute the PSD using a fast Fourier 
transform algorithm as points are equidistant in time in the light 
curve. The result is plotted in Fig. [3] (grey line). The PSD re- 
binned over 12.8 //Hz (101 bins) and over the large separation 
are also shown (black and green lines respectively in Fig. [3j. 

A single p-mode bump is clearly visible corresponding to 
HD 169392A. The convective background (red curve in Fig. [3]l 
can be fitted using a standard maximum likelihood estimator 
(MLE) with 3 standard components: a white noise (W), one 
Harvey law (|Harvey 1985), and one power law (see for more 
details, e.g., |Mathur et al.|20lTb) : 



B(v) = W + 



4j g a\ 



1 + (27TT„v) ff s 



+ P„v- 



(1) 



where t ? and cr g are the characteristic time scale and amplitude 
of the granulation, and a g and e a are the exponents character- 
izing the temporal coherence of the phenomenon. The obtained 
coefficients are: W = 0.945 + 0.003 ppm 2 //^Hz, r g = 557 ± 60 
s, o- g = 46.8 + 2.1 ppm, a g = 2.43 + 0.18, P a = 24 + 4.3, and 
e a = 1.8 ± 0.06. The value obtained for the granulation time 
scale agrees with the relation derived with the Kepler red giants 
with v max (Ma thur et al.|20lTb| . Finally, we also computed the 
rms amplitude of the granulation: 



(2) 



and we obtain cr gl - ms ^43.3 ±2.2ppm. Nevertheless, we have 
to keep in mind that this measurement is polluted by the sec- 
ondary star of the binary system. There are two opposite effects. 
HD 169392B contributes for 20% to the total observed flux (see 
also Sect. 7.6). This translates into an underestimation of cr grms 
by a similar amount. On the other hand, the granulation com- 
ponent measured for HD 169392A probably includes a contri- 
bution from the granulation (and other short-frequency variabil- 
ities) of HD 169392B. This would lead to an overestimation of 



4.2. Looking for HD 169392B 



5.2. Surface rotation 



We derived in Section |2.1| that the maximum of the p-mode 
bump for HD 169392B should lie around a region centered at 
~ 2650 //Hz. Five different teams looked for the signature of this 
second star in a blind way, without success. We failed to obtain 
a possible detection at a level of 90 % confidence level in all the 
cases. With the stellar parameters of the two components {T^, 
my, R), we computed the probability of detection of the modes 
following the method described in Chaplin et al. ( |2011a I. For the 
A component, we obtain a 100% probability of detection. For 
the B component, which has its flux dimmed by the A compo- 



Magnetic features such as starspots crossing the visible stellar 
disk of a star produce a fluctuation of the flux emitted by the 
star. The careful analysis of these long-period modulations in 
the observed light curve provides invaluable information on the 
average rotation rate of the surface of the star, at the latitudes 
where these magnetic features evolve. This study can be done 
directly in the light curve by modeling the spots (e.g. Mosser 



et al. [2009b a ), or by the close inspection of the low-frequency 



part of the power spectrum (e.g. |Garcia et al.||2009 Campante 
etal.|2011 " 
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Fig. 3. PSD of HD 169392 data (grey) that has been modeled as 
explained in the text (red curves) using the usual three compo- 
nents: white noise (W), granulation noise (B g ) and stellar activity 
and/or large scales of convection (B a ). The red continuous line 
is the sum of the three components above mentioned, resulting 
from the fit of the smoothed spectrum over Av (green line). The 
black curve corresponds to the spectrum smoothed over 101 bin, 
i.e. 12.8juHz. 

Looking at the light curve displayed in Fig. [2] we do not 
see any clear modulation produced by spots. However, the low- 
frequency part of the power spectrum is dominated by two sig- 
nificant peaks at 0.25 and 0.64 pMz (see black line of Fig. [4}, 
the latter being a harmonics of the first peak when we con- 
sider the frequency resolution of 0.127 pHz. But as explained 
in Section[3] there are two jumps of possible instrumental origin 
at 43.3 and 80.05 days that could be at the origin of the peak 
at 0.25 yuHz. To verify this hypothesis, we correct the jumps us- 
ing the same procedure used to correct Kepler data ( |Garcia et al. 
|2011) . In the PSD of the resultant light curve (blue line inFig.|4J 
there is no power at the above mentioned frequencies. Therefore, 
it is probable that the 0.25 and 0.64 yuHz peaks have been pro- 
duced by instabilities in CoRoT but a stellar origin could still be 
possible. 

6. Methodology followed to extract the mode 
parameters of HD 169392 A 

Nine teams estimated the mode parameters for HD 169392 A 
with two teams providing two sets of results from differing meth- 
ods. This gives a total of eleven sets of mode parameters. We 
describe, albeit briefly, the methods used by the different teams 
to fit the power spectrum corrected from the background. A list 
of the methods used for each of the sets of mode parameters can 
be found in Table [4] 

All teams used a common philosophy for estimating mode 
parameters: maximise the likelihood function for a given model. 
The models used by the teams are common in that all are com- 
posed from a sum of Lorentzian profiles to describe the modes 
of oscillation, plus a background to account for the instrumental 
and stellar noise (see Section [5T| . Each team used the frequency, 
height, linewidth, rotational splitting, and angle of inclination to 
describe a mode of oscillation (e.g. |Appourchaux et aL"1 |2008). 
All the teams fitted a single linewidth for each order. The main 
differences between teams arises from the method used to max- 
imize the likelihood. The three broad categories used are max- 




2 3 
Frequency (/xHz) 

Fig. 4. Low-frequency part of the power spectrum between 0. 1 
and 5 pMz. The black line correspond to the light curve showed 
in Fig [2] Two peaks are visible at 0.25 and 0.64 yuHz. The blue 
line correspond to the PSD of the light curve in which the jumps 
at 43.3 and 80.05 days have been corrected. 

Table 4. Description of the fitting methods of each team. 



Fitter ID 



Method 



Splittings Angle 



OB° 
JB b 
GRD* 
RAG* 
RAG* 
SIT 
RH rf 
CR* 
CR* 
DS e 
GAV* 



MCMC global Free Free 

MAP global Free Free 

MAP global Free Free 

MLE global Free Free 

MLE global Fixed Fixed 
MLE pseudo-global Free Free 

MCMC global Free Free 

MLE global Free Free 

MLE global Fixed Fixed 
MLE local Free Free 
MLE global Free Free 



(fl) |BenomarN2008t 
m |Appourchaux et al.|<fl998} 
^Fletcher etal.N201 It 
w |Handberg & CampantelfcOll) 
w |SalabertetaL| ( |2004{ 



imum likelihood estimation (MLE, Appo urchaux et ah] 1998| l, 
that has a variant called the maximum a posteriori method (MAP, 
I Gaulme et al.|2009|>, and Markov Chain Mon te Carlo (MCMC, 
|Benomar|2008||Handberg & Campante|201 1| >. 
The MLE or the MAP approaches were used by the majority 
of teams. Briefly, MLE estimates the mode parameters by per- 
forming a multi-dimensional maximisation of the likelihood pa- 
rameter. Typically formal one-sigma errors are determined from 
the Hessian matrix. Like all Bayesian approaches, MAP tech- 
niques extend this concept by including a priori knowledge by 
implementing penalties on the likelihood parameter for defined 
regions of parameter space. As for MLE, MAP errors are taken 
from the Hessian matrix. Both MLE and MAP can be applied 
either globally, that is simultaneously fit to the entire frequency 
range (e.g. Appourchaux 2008), or locally, considering parame- 
ters separately over one large separation ( Salabert et al.| [2004). 
When coupled with MCMC, a Bayesian approach is a power- 
ful way (but also time consuming) to extract the full statistical 
information with respect to the likelihood and our prior knowl- 
edge by mapping the probability density function (PDF) of each 
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parameter. Contrary to MLE or MAP, MCMC techniques avoid 
to be trapped in local maxima and no assumptions on the na- 
ture of the probability density function are needed. Thus such 
approach provides more conservative and robust results, espe- 
cially for the uncertainties as their estimation do not rely on 
the Hessian matrix, but on the cumulative distribution function, 
built thanks to the PDF. In the present paper, two fitters (OB, 
RH) used MCMC approaches with different priors. However, a 
careful choice of the priors was needed as it may have a strong 
impact on the best fitted parameters, especially at low signal- 
to-noise ratio. This star is a perfect example: £ — and £ = 2 
frequencies lower than 870 /iHz (the three first radial orders) 
were barely fitted (posterior probability highly multimodal and 
non Gaussian) if no smoothness condition was applied onto the 
frequencies of these modes. A smoothness condition acts as a 
filter and avoid strong, erratic variations of frequencies from one 
order to another (Benom ar et al.|20T2) . 

Each of the three approaches can be applied in slightly differ- 
ent ways. Two fitters (RAG & CRR) produced two (MLE/MAP) 
results with and without the two parameters describing rota- 
tion rate and angle of inclination. In all cases for this analysis 
mode profile asymmetry was not included. This is due to insuffi- 
cient signal-to-noise ratio: asymmetry fitting requires very good 
signal-to-noise ratio and high resolution. 

The analysis has revealed the presence of mixed modes. As a 
consequence, to complement the fitting techniques, we have also 
used the asymptotic relation found by Goupil (private communi- 
cation), following the ideas originally developed by Unno et al. 
1989) and tested by Mosser et al. ( 2012| l. This methodology is 



particularly suited to identify the presence of mixed modes in 
the power spectrum and to derive the measure of the gravity pe- 
riod spacing. This global seismic parameter provides a strong 
constraint on the core radiative region of the interior models. 

Given eleven sets of fitters results, we required a statistically 
robust method of comparison to produce a final list of mode 
parameters for future modelling. Here we adapted a previously 
used method (Appourchaux et al. 2012c) to determine a con- 
fident list of mode frequencies. This method used an iterative 
outlier rejection algorithm to populate two lists (a maximal and 
a minimal list) of modes which are then used to determine the 
"best" fitter who will provide the frequencies and is defined as 
the fitter with the smallest normalized rms deviation from the 
average frequency. 

At each radial order, n, and degree, /, we compared N different 
estimated frequencies. As a first cut a visual inspection of all 
mode parameters and errors was performed. Modes with clearly 
incorrect parameters (for example, very small/large linewidths or 
very large errors) were rejected. Statistical outlier rejection was 
performed with Peirce's criterion ( Pei rce| 1 852| [Gould 1855), a 
method that is based on rigorous arguments. Peirce's criterion 
was applied iteratively to each mode frequency set until no data 
points are rejected. The method can be described with the fol- 
lowing pseudo-code: 

- COMPUTE mean x and rms cr deviation from the sample x. 

- COMPUTE rejection factor r from |Gould| ( |1855[ l assuming 
one doubtful observation. 

- REJECT: Reject data if \x t -x\> rcr. 

- IF n data are rejected THEN compute new r assuming n + 1 
doubtful observations ELSE END. 

- GOTO REJECT. 

With the results of the outlier rejection we populated the max- 
imal and minimal list in the following way. The minimal list 
was designed such that we had a very high level of confidence in 



the results included. Previous works have used varying threshold 
criteria for inclusion in the minimal list. Metca lfe et al.| ( |2010| l 
required that all fitters agree to place a frequency in the mini- 
mal list, i.e. N results are accepted. In an updated version of the 
recipe, ( Campante et al. 201 1] > and Mathur et al. ( |201 la} required 
that N/2 results are accepted. Here we followed the N/2 accep- 
tance criterion for the minimal list. Modes forming the maximal 
list must have at least two results accepted by Peirce's criterion, 
because of this the maximal list must be treated with some cau- 
tion. 

In the maximal and minimal lists we selected the actual frequen- 
cies of the modes in each list. We calculated the normalized root- 
mean-square (RMS) deviation of the N teams with respect to 
the mean of the frequencies for each mode in the minimal list 
(see Equation 2 of Mat hur et al.||201 la) . The data set with the 
lowest normalized RMS deviation was declared the "best" fitter. 
Two methods (MLE and MCMC) obtained very similar values 
of RMS and we decided to present in the following sections the 
results obtained by the MCMC, providing a reproducible set of 
frequency parameters. 



7. Results 

Following the methodology described in the previous section, 
we obtained the mode parameters presented in this section. 



7. 1 . Mode frequencies 

The comparison of the 11 sets of the modes frequencies pro- 
vided by the 9 teams gave the minimal and maximal lists of fre- 
quencies. The fitter selected from the smallest RMS deviation 
(namely OB) used the MCMC algorithm and the model fitting 
the data is represented on Fig. [5] We obtained a minimal list of 
25 modes. The maximal list contains 4 additional modes. The 
lower panel of Fig. [5] shows the ratio of the power spectrum to 
the fitted model. The power spectrum can be thought as com- 
posed of a limit spectrum multiplied by a x 1 with a 2 degree- 
of-freedom noise function (e.g. Anderson et al.|l 990) and when 
smoothed over k bins the noise function is distributed as a x 1 
with 2k degrees of freedom (e.g. Appourchaux et al. 1998 ). If the 
fitted model describes the limit spectrum well, then the ratio of 
the power spectrum to the model will be distributed as the noise 
function. Hence, the power spectrum normalized by the model 
in the lower panel of Fig. [5]reveals structures in the power spec- 
trum not compensated by the fitted model in the region where 
the £ — 3 modes are expected (shown with arrows in the plot). 



7.2. 1=3 modes 

With the MCMC technique, we fitted the £-0, 1, 2, and 3 si- 
multaneously. We checked that the frequencies of the modes £ 
=0, 1, and 2 remain unchanged. In Table [5] we give the minimal 
and maximal lists of frequencies as well as the three £ = 3 modes 
that have a significant posterior probability. However, we remind 
the reader that these modes have to be taken very cautiously. By 
adding the fit of the I — 3 modes, it slightly modified the final 
values of the amplitudes, widths, splittings, and inclination an- 
gle. The values listed in Table [7] correspond to the fitting of the 
four ridges. 
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Fig. 5. Upper panel: Power spectral density (PSD) of HD 169392 in the p-mode region at full resolution (grey) and smoothed over 
15-bin wide boxcar (black). The red line corresponds to the fitted spectrum from refit after the frequency comparison stage. Lower 
panel: Power spectral density for the p-mode region divided by the fitted model at full resolution (grey) and smoothed over 15 
bins (black). Regions highlighted in green cover the regions corresponding to the harmonics of the CoRoT orbital period such that 
v - (/ix 161.7) ± (m X 1 1.5) //Hz, where n and m are integers. The arrows show the positions where we would expect ( = 3 modes. 



7.3. Mixed modes 

By looking at the echelle diagram (Fig.[6]i, we notice an avoided 
crossing, corresponding to the presence of a mixed mode (e.g. 
|Osaki|1975||Aizenman et al.|1977HDeheuvels & Michel |2011| i. 
The frequency pattern is affected by the coupling of the pressure 
wave in the outer convective region with a gravity wave in the 
radiative interior. Hence, we can use the formalism developed 
for mixed modes in giants (Mos ser et al.|2012] l to describe this 
pattern. 

For p modes, the second order correction of the asymptotic 
relation is related by a curvature term a: 



K + 2 



£ P - doe + ^ K 



] 2 Av, 



(3) 



where e p is the pressure phase offset, doc accounts for the small 
separation, and n max = v max /Av. The g modes are reduced to the 
first order term: 



+ £ g ) Ani 



(4) 



where n g is the gravity radial order, e g the gravity phase offset 
and Ani the gravity period spacing of dipole modes. 

The coupling of the p and g waves is then expressed by a sin- 
gle constant (Moss eret al.|2012 i. Table [6] shows the asymptotic 
parameters used for the fit. Compared to modes in giants that 
have very large gravity orders (in absolute value), mixed modes 
in a subgiant with Av ^ 56/vHz have small |n g |, Therefore, we 
cannot consider the asymptotic constant e g to be 0. 



Table 6. Mixed-mode parameters 



Asymptotic p modes 



Pressure phase offset 
Small separation 
Curvature 



1.17 ±0.02 
0.03 ±0.02 
0.0029 ± 0.0008 



Asymptotic g modes 



Gravity period spacing 
Gravity phase offset 



AIL 



476.9 ± 

0.07: 



4.3 s 

:0.03 



Mixed modes 



Coupling 



0.098 ±0.005 



The best fit of the radial ridge, obtained from a global 
agreement over 9 orders, is provided with Av = 56.3 yuHz and 
s p = 1.17. We then infer a gravity period spacing of the or- 
der of 476.9 + 4.3 s, with a precision limited by the uncertainty 
on the unknown gravity offset. Best fits are obtained with e g = 
0.07 + 0.03. In all cases, the coupling is small: q = 0.098 ±0.005. 
Two mixed modes are reported at 816 and 1336 ^Hz with grav- 
ity orders of -2 and - 1 , respectively, corresponding to the period 
of g modes (the avoided crossing), P g = QnJ + 1/2 + e g ) An! 
(Goupil et al. in prep.). These mixed modes bring strong con- 
straints on the age of a star as the coupling between the p- and 
g-mode cavities varies very fast with the stellar age, introduc- 



ing strong constraints on the age of the star (e.g. 



Christensen-| 



Dalsgaard |2004l |Metcalfe et al.|2010[ |Mathur et al.|2012) . The 
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10 20 30 40 

Frequency modulo 56.30 //Hz 

Fig. 6. Echelle diagram of the CoRoT target, HD 169392A, with 
the minimal (close symbols) and maximal (open and filled pink 
symbols) lists of frequencies, and the frequencies from the best 
fit model of AMP (white symbols): I = (circles), I = 1 (trian- 
gles), t = 2 (squares), and I — 3 (stars) as in Section 7.2. The pink 
triangles without uncertainties correspond to the mixed modes 
(^=1) mentioned in Section 7.3. 

presence of the mixed modes in HD 169392 A indicates that this 
star is quite evolved. 

Unfortunately, none of these two modes were fitted by the 
9 teams. The peak at 816//Hz is surrounded by the orbital har- 
monics and without a proper modeling of these peaks, it is not 
possible to fit the mixed mode. The high frequency mixed mode 
at 1336 /iHz has a too weak signal to be fitted. 

7.4. Frequency differences 

We fitted the frequencies of the t — modes against the order n 
and according to Eq. 3, the slope is the mean large frequency sep- 
aration. We obtained (Av)=55.98 + 0.05 yuHz, which agrees with 
the value obtained from the global techniques (see Section|4| and 
the mixed-mode asymptotic fit (Section 7.3). The variation of the 
mean large separation is shown in Fig. [JJ We notice the dip for 
the modes 1-1, due to the avoided crossing already mentioned. 
The small oscillation we observe is typical of the signature we 
can have for the helium's second ionization zone below the stel- 
lar surface. We would need more modes to be able to extract the 
position the second ionisation zone. 

We also computed the variation of the small separation be- 
tween t — and 2 following Roxburgh & Vorontsov ( 2003 i and 
it is shown in Figure [8] The small separation presents a general 
decreasing trend (with a slope of -0.007), which is very simi- 
lar to what we observe on the Sun. The mean value is around 
4.14/jHz. This value fits perfectly the CD diagram presented 
by |White et al.| ( [201 l| l where for subgiants, the small separation 
reaches approximately similar values for stars with very different 



7.5. Linewidths 

The linewidths of the modes (T( A ) are related to their lifetimes 
(Tt,n) following the formula for I — and order n: 
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Fig. 7. Large separation as a function of frequency for t — 0, 
1 and 2. The dashed line represents the mean large separation 
of the radial modes. The red squares at low and high frequency 
represent the values obtained using the modes of the maximal 
list. 
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Fig. 8. Small separation, <Jo2> as a function of frequency. The red 
squares at low and high frequency represent the values obtained 
using the modes of the maximal list. 



The results from the MCMC method, which fitted a single 
linewidth for 1, and 2 modes per order, are given in Table [JJ 
and represented in Fig. [9] We note that the trend is rather flat with 
a median value of 1.8 //Hz, with a little increase with frequency. 
The first and last points were obtained with the frequency of the 
maximal list of modes. This behavior is similar to the one found 
in the Sun in the plateau region and the higher frequencies. We 
also observe a small oscillation followed by a small increase of 
the linewidth, which is also seen in other G-type stars observed 
by CoRoT (e.g. |Deheuvels et al.|2010"}|Ballot et al.|2011b| > . 

Several works have been done to show the dependence of the 



linewidths with the effective temperature of the star as F ~ r eff . 
Cha plin etal.| ( 2009) derived a relation where the mean linewidth 



scales a s for temp eratures ranging between 6800 K and 



T(),„ = (7rFo ; „) 



(5) 



5300 K. Baud in et al.| ( |2011| l re-evaluated the exponent to be 
16 + 2 using 5 main-sequence stars observed by CoRoT and 
the Sun. The latest work by Appourchaux et al. ( 2012a| l based 
on the analysis of 42 cool main-sequence stars and subgiants 
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Fig. 9. Linewidths of I = modes as a function of frequency. 
Same legend as in Fig. [8] A single linewidth has been fitted for 
each order. 



observed by Kepler refined the expression to F ~ T!Z . For 
HD 169392, we measure the linewidth at maximum height as 
1.41 //Hz (€ = 0, n = 18) that agrees very well with Fig. 2 of 
App ourchaux et af] (2012a) using T e g = 5885 K and with the the- 



oretical computation by |Belkacem et al. (2012 1. Finally, [Corsaro 
et al. (2012 1 analysed the linewidth of red-giant stars in clus- 
ters and found that by taking into account main-sequence and 
subgiant stars, F varies exponentially with T e ff), which gives 
1 .66 //Hz, still in agreement with the observed value. 



7.6. Amplitudes 

From the heights and linewidths of the modes, we can compute 
the rms amplitude according to the following expression: 



(6) 



where T^„ is the width of the mode (€, n) and H( n is the 
height of that mode in a single-sided power spectrum (hence 
the factor 2 in the numerator). The values obtained are listed 
in Table [7] and shown in Fig. 10 We can see the profile of 
the p-mode bump. The maximum amplitude of the modes is of 
3.82 +0.30 ppm, which agrees with the value obtained in Section 
4.1 (3.61 +0.35 ppm). Due to the imperfect observational win- 
dow (with a duty cycle of ~ 0.83), there is a non-negligible leak- 
age of mode power in small aliases. We need to correct the 
amplitude for this effect and this leads to an observed value 
of 4. 12 + 0.32 ppm. To take into account the response function 
of CoRoT, we convert the amplitude to the bolometric one by 
computing /? osc =7.87 and R e=0 =4A4 (see |Michel et al.|2009] for 
details), which gives the bolometric correction Cb o i=0.90 ± 0.01. 
Finally, since HD 169392 is a binary, a part of the observed total 
flux comes from the secondary component. With the magnitude 
quoted in Section 1, we find L A /L B =3.9, which allows us to 
conclude that ~ 20% of the flux is dimmed by the second com- 
ponent. Therefore, the final value for the maximal amplitude of 
the modes is A^ x = 4.45 + 0.34 ppm. 

Several expressions have been developed the last few years 
to compute the maximum amplitude of the modes. According 



Fig. 10. Amplitude of £ = modes as a function of frequency. 
Same legend as in Fig. [8] 



to Samadi et al. (2007b), the theoretical value of the maximum 
amplitude is obtained with: 



/ L/Le 



\M/M a ) \ T t 



V 5777 

I \ ~z. A c . 



(7) 



with 5=0.7 and according to Michel et al. (2009), 
A0,max=2.53 + 0.1 1 ppm (in rms). Using the values of R 
and M from the AMP model (see Section 8), we obtain 
6.02 ±0.60 ppm. 

We computed the theoretical value with the updated scaling 
relation derived by Kjeldsen & Bedding (201 1 ), which depends 
on the lifetime of the modes, t osc : 



u 



.0.5 



{M/M o y- 5 T 



1.57-2.25+r 



(8) 



eff 



We took r=1.5 as adopted by Michel et al. (2008! and we 



obtain 4.32 ± 1.10 ppm, which is much closer to the observed 
amplitude. 

Huber et al. ( 201 1[ ) used several hundreds of stars (red gi- 
ants and main-sequence stars) to derive an empirical relation, 
slightly different from the previous one and that is similar to the 
results obtained from red giants observed in clusters ( Stello et al. 
201 la) . This relation is more dependent on the mass of the star: 



(L/L Q y 5777 

(M/M Q )' T eS 1 



(9) 



with 5=0.838, f=1.32, and A 0max =2.53 ±0.11 ppm (in rms), 
leading to 6.58 ±0.79 ppm for|Huber etaL] (|20 1 1 1>. |Stello et al. 



( 201 la| l found slightly different parameters: s=0.95 and f=l. 
(in the adiabatic case), which gives A max = 7.21 ± 0.92 ppm. 

The observed value is at the lower limit of the theoretical 
ones within 3-cr for modelings in equations (7) and (9) and 
agrees within 1-cr for the prediction from equation (8). 

7.7. Relative visibilities 

Accurate estimations of the relative visibilities of non-radial 
modes are important for our knowledge of the physics of the 
stellar atmospheres as they depend on the stellar limb darken- 
ing and the observed wavelengths. Although well defined in the 
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Sun at different heights in the solar photosphere (Salabert et al. 
|201 1) , their measurements are more critical for other stars due to 
lower signal-to-noise ratio and shorter time series. Nevertheless, 
the mode visibilities were successfully measured in two of the 
CoRoT stars: HD 522 65 ( |Ballot et al.||2011b| l and HD 49385 
(Deheuvels et al. 20 1 0] > . The uncovered visibilities quantitatively 
agree well with the theoretically calculations of limb-darkening 
functions from Ballot et al. (201 la i, who also showed a depen- 
dence on T e ff . 

Estimations of the relative visibilities, Vf, of(- 1,2, and 3 
modes were obtained for HD 169392 A. Consistent values within 
the error bars were obtained among the different peak-fitting 
techniques. Some tests using different inclination angles and 
rotational splittings showed that the uncovered visibilities give 
similar results within the uncertainties. The relative visibility of 
the modes I — 1,2, and 3 measured by the MCMC analysis are 
given in Table [8] as well as the median values of their distribu- 
tions and their associated l-cr values. 



7.8. Splitting and inclination angle 

We took the step to check that mode parameters were indepen- 
dent of the assumed angle of inclination. Multiple fits with a 
fixed rotational splitting range of fixed angles were tested based 
on the MCMC results. The mode frequencies returned by the 
various fits were virtually independent of the angle used. A 
last check was done by directly looking at the correlation map 
between width, inclination angle and splittings from the MCMC 
technique. The mode linewidths and amplitudes returned by 
these fits showed a small dependence on the angle of inclination 
and the splittings, but at a much lower level than the l-cr error 
bars over the range 15 ° < i < 85 ° - the range limits suggested 
by MCMC results. 




40 60 
Inclination (°) 

Fig. 11. Correlation map between rotational splitting and incli- 
nation. The red (resp. yellow) region correspond to the l-cr (resp. 
2-cr) confidence area. 

Table 8. Medians and standard deviations for the inclination 
angle (i), the rotational splitting (v s ), the projected splitting 
(v s sin i) and the relative visibilities of the £ — 1,2, and 3 modes 
(Vi=i, Vi=2, V/=3) obtained from MCMC. These statistical infor- 
mation are deduced from the joint-probability density function 
of the inclination/splitting and from the probability density func- 
tion for the other parameters. 
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Fig. 1 1 shows the correlation map of the splitting and the in- 
clination (two-dimensional probability density function). We no- 
tice two maxima of almost equal heights into the joint-posterior 
probability density function. This bi-modality suggests that two 
different values for the joint-parameters are possible. However, 
the two solutions are compatible within 2 cr and therefore are not 
statistically different (cf. Table [8}. In order to have a more strin- 
gent determination, one would need additional constraints, such 
as the signature of the surface rotation, unfortunately unavailable 
with the present dataset (see Section 5.2 1. 

The product of the two quantities, which gives the projected 
splitting, v s sin i, is well defined as seen by the curved ridge and 
the uncertainties presented in Table [8] However, as anticipated 



2008), which is also used to build the standard model of the 



Sun. It includes the OPAL 2005 equation of state (Rogers & 
Nayfonov 2002) and the most recent OPAL opacities (Iglesias 



& Rogers p996|l associated with the low-temperature opacity 



table of Alexander & Ferguson (1994). Finally, the convection 



is treated according to the mixing-length theory (Bohm-Vitense 



by |Ballot et aL] (2006 



the robust estimation of the individual 
quantities of angle and splitting is hampered by two competing 
solutions, the first with relatively low angle and high splitting, 
and the second with high angle and lower splitting. Our inability 
to discern between the two solutions prevents us from reporting 
the desired 'robust' result. 



8. Modeling HD 169392A 

With the spectroscopic and asteroseismic information retrieved 
in the previous sections, we modeled HD 169392A with the 

2009). 



Asteroseismic Modeling Portal (AMP, Metcalfe et al. 
AMP is a genetic algorithm that aims at optimizing the match be- 
tween the observed and the modeled frequencies of the star in ad- 
dition to the spectroscopic constraints. The models are computed 
by the Aarhus STellar Evolution Code (Christensen-Dalsgaard 



1958), and diffusion and gravitational settling of helium follow 
the prescription by Michaud & Proffitt] ( 1993 i. A surface cor- 
rection was applied to the frequencies following the prescription 
by |Kjeldsen et al.| ( |2008b]>. AMP has already modeled a large 
number of stars (e.g . |Metcalfe et al.||2010[|Mathur et al.||2012 
|Metcalfe et al.|2012] > for a given physics. 

For the observables, we used the spectroscopic constraints 
from Table [T] and the mode frequencies from the minimal list 
(Table [51. The best fit model obtained has a x 2 of 102. Since 
seismic and spectroscopic observables are used, we also com- 
puted two separate normalized^ 2 : ^ 2 pec =0.22 and ^ 2 c(j = 1.43. 

The xlpec shows a good agreement between the observables 
and the best fit model. In Figure [6] the modeled frequencies are 
superimposed on the echelle diagram obtained with CoRoT, and 
we can see that the best fit model frequencies match well the ob- 
servations. We can also notice that the maximal list frequencies, 
the I — 3 modes included, are well reproduced by the model, sug- 
gesting that these frequencies probably correspond to real signal. 

We obtained a mass of M = 1.15 + 0.01 Mq, a radius of 
R= 1.88 ± 0.02 Rq, and an age t = 4.33 ±0.12 Gyr, where the 
uncertainties quoted are internal errors. 
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Fig. 12. Evolutionary track computed by AMP and position of 
the best fit model with the spectroscopic uncertainty (blue sym- 
bol). 



Using the parallax n = 13.82 ±1.2 mas, we obtained a lumi- 
nosity L = 4. 16 ± 0.6 L G while the value from the best-fit model 
of AMP is 4.18 ± 0.01 L , which is a very good agreement. 
We have also followed the methodology developed by |Silva 



Agui rre~eTaLl ( pOT2" 



to couple the Casagrande et al. 2010 im- 
plementation of the InfraRed Flux Method (IRFM) with astero- 
seismic analysis. These technique provides results using both the 
direct and grid-based methods, as well as a determination of the 
distance to the target. The results are in complete agreement with 
the AMP ones for the loggand T e ff (see Table |5J. The mass and 
radius agree with the first method within 1 <x and with the sec- 
ond one within 2 cr. The distance determined by these techniques 
is d=70± 5.3 pc, in excellent agreement with parallax measure- 
ments. 

The best fit model of AMP confirms that this is indeed a 
subgiant that has exhausted its central hydrogen and has no con- 
vective core. As seen on the HR diagram with the evolutionary 
track for this model (Figure 12 1, it has not yet started ascending 
the red giant branch. 



9. Conclusion 

In the present analysis, we have derived the fundamental param- 
eters of the two components of the binary system using HARPS 
observations. We obtained very close values of metallicity for 
the two components. This is quite similar to what has been 
observed for other binary systems (e.g. 16Cyg Metcalfe et al. 
2012| l, suggesting that the A and B components were formed 
from the same cloud of gas and dust. None of the two compo- 
nents have significant amount of Li. 

The analysis of 91 days of the CoRoT observations did not 
allow us to detect a significant seismic signature of the B com- 
ponent. This result agrees with the low detectability probability 
of 2% for the B component. 

We measured the global parameters of the modes of 
HD 169392A leading to <Av> of 56.98 ± 0.05 fMz and v max of 
1030±55jt/Hz making this star very similar to another CoRoT 
target, HD 49385 ( |Deheuvels et al.|2010) . 

Nine groups reported p-mode parameters and the compari- 
son of the results led to a minimal and maximal lists of frequen- 
cies. Some significant power was found in the region where we 
would expect I — 3 modes. We fitted the t — 0, 1,2, and 3 modes 



all together with the MCMC technique and obtained that three 
modes t — 3 are significant according to the posterior probability. 
We analysed the mixed modes of the star using their asymptotic 
properties and detected two avoided crossings at 816//Hz and 
1336/iHz. From these avoided crossings, we derive the grav- 
ity spacing AHi =477 + 5 s. We conclude that this star is quite 
evolved but not yet on the red-giant branch. 

We derived the amplitude of the modes, which is smaller 
than the predicted amplitudes within 1- to 3-cr depending on 
the theoretical formula applied. The linewidths of the modes ob- 
tained are coherent with the dependence on the effective temper- 



ature derived by |Appourchaux et al.| ( |2012a| ) and |Corsaro et al 
( |20T2) . 

The study of the splittings and inclination angle led to two 
possible solutions. Either we obtain an inclination angle of 
20 - 40 ° with splittings of 0.4 -1.0 //Hz, or we have an inclina- 
tion angle of 55 - 85 ° and splittings of 0.2 - 0.5 //Hz. We note 
that no clear signature of the surface rotation has been detected 
and in absence of additional constraint, the significant correla- 
tion between the splitting and the inclination angle prevents us 
from constraining them more precisely. 

The global parameter of the modes allowed us to have a first 
estimation of the mass and radius of the star using scaling rela- 
tions based on solar values: 1.34 + 0.26 M G and 1.97 ± 0.19 Rq. 
A more thorough modeling with AMP, based on the match of 
the individual mode frequencies and the spectroscopic parame- 
ters obtained in this work, provided M = 1.15 ± 0.01 M , a radius 
of R = 1.88 + 0.02 R , and an age r = 4.33 ± 0. 12 Gyr. 
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Table 5. Minimal and maximal lists of frequencies for HD 169392 in fiHz obtained from 91 days of CoRoT observations with 
MCMC. 



Order 


1 = 





t=\ 


e = 2 


e = 3 


12 








743.22 ± 1.34° 




13 


748.35 ± 


0.40" 


770.96 ±0.37 


797.97 ±0.76 




14 


804.50: 


t0.35 


838.21 ±0.38 


853.50 ±0.85 




15 


859.70: 


t0.38 


886.83 ±0.23 


909.34 ±0.46 




16 


914.14: 


t0.27 


940.52 ±0.27 


965.45 ± 0.41 


930.52 ± 1.87* 


17 


969.77 : 


t0.46 


996.40 ±0.32 


1022.18 ±0.27 


986.08 ±1.20* 


18 


1026.79 


+ 0.27 


1052.68 ±0.20 


1080.04 ±0.44 


1043.43 ± 1.23* 


19 


1083.77 


±0.39 


1109. 19 ±0.36 


1137.05 ±0.58 




20 


1139.76 


±0.42 


1166.17 ±0.34 


1193.29 ±0.97 




21 


1196.71 


±0.73 


1221.76 ± 1.38 


1248.79 ± 1.44" 




22 


1253.73: 


tO.75" 








{a) Modes that belong to the maximal list. 



Significant modes according to the MCMC method. 



Table 7. Amplitudes and linewidths for all the modes of Table plcomputed with the MCMC. 



/ 


n 


Linewidth (/jHz) 


+<T 


—cr 


Amplitude (ppm) 


+cr 


—cr 





13 


0.80 


1.02 


0.64 


1.63 


0.33 


0.32 





14 


1.39 


0.65 


0.62 


2.40 


0.30 


0.32 





15 


1.35 


0.39 


0.30 


2.83 


0.29 


0.27 





16 


1.30 


0.63 


0.43 


2.81 


0.31 


0.26 





17 


2.37 


0.61 


0.52 


3.82 


0.30 


0.29 





18 


1.00 


0.57 


0.27 


3.72 


0.36 


0.37 





19 


2.01 


0.77 


0.54 


3.39 


0.31 


0.29 





20 


1.97 


0.58 


0.47 


3.23 


0.24 


0.25 





21 


4.14 


1.53 


1.18 


2.54 


0.25 


0.27 





22 


1.27 


1.35 


0.77 


1.36 


0.29 


0.36 




13 


0.80 


1.02 


0.64 


1.94 


0.38 


0.40 




14 


1.39 


0.65 


0.62 


2.84 


0.37 


0.37 




15 


1.35 


0.39 


0.30 


3.37 


0.32 


0.31 




16 


1.30 


0.63 


0.43 


3.35 


0.34 


0.32 




17 


2.37 


0.61 


0.52 


4.53 


0.34 


0.33 




18 


1.00 


0.57 


0.27 


4.40 


0.42 


0.38 




19 


2.01 


0.77 


0.54 


4.04 


0.33 


0.33 




20 


1.97 


0.58 


0.47 


3.83 


0.31 


0.31 




21 


4.14 


1.53 


1.18 


3.02 


0.31 


0.33 


2 


12 


0.80 


1.02 


0.64 


1.34 


0.28 


0.27 


2 


13 


1.39 


0.65 


0.62 


1.96 


0.28 


0.27 


2 


14 


1.35 


0.39 


0.30 


2.33 


0.25 


0.24 


2 


15 


1.30 


0.63 


0.43 


2.32 


0.25 


0.23 


2 


16 


2.37 


0.61 


0.52 


3.13 


0.25 


0.25 


2 


17 


1.00 


0.57 


0.27 


3.04 


0.32 


0.29 


2 


18 


2.01 


0.77 


0.54 


2.79 


0.26 


0.26 


2 


19 


1.97 


0.58 


0.47 


2.64 


0.27 


0.24 


2 


20 


4.14 


1.53 


1.18 


2.08 


0.23 


0.25 


2 


22 


1.69 


1.42 


0.91 


0.99 


0.27 


0.29 


3 


16 


1.30 


0.63 


0.43 


0.65 


0.07 


0.06 


3 


17 


2.37 


0.61 


0.52 


0.89 


0.07 


0.07 


3 


18 


1.00 


0.57 


0.27 


0.86 


0.08 


0.09 



Table 9. Modeling results from AMP, direct method + IRFM, and grid-based method +IRFM. The uncertainties for the AMP best 
fit model are the internal errors. 



Method 


M(M G ) 


R (R Q ) r (Gyr) 


logg(dex) 


T eB (K) 


L/L Q 


AMP 


1.15 + 0.01 


1.88 ±0.02 4.33 ±0.12 


3.95 ±0.01 


6016 ±5 


4.18±0.01 


Direct method + IRFM 


1.30 ±0.23 


1.93 ±0.13 


3.97 ±0.02 


6002 ± 77 


4.45 ± 0.67 


Grid-based method + IRFM 


i it +0.05 
-0.03 


1.96 ±0.02 


3.95 ±0.01 


6002 ± 77 


4.51 ±1.04 
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